The Architecture of Higher-Order ODEs
An $n$th-order linear differential equation is characterized by its highest derivative. We define the general form as Equation (1):
$$P_0(t) \frac{d^n y}{dt^n} + P_1(t) \frac{d^{n-1} y}{dt^{n-1}} + \dots + P_{n-1}(t) \frac{dy}{dt} + P_n(t)y = G(t)$$ (1)
To facilitate theoretical analysis, we often normalize this equation by dividing by $P_0(t)$, assuming it is non-zero on the interval of interest. This yields the Standard Form (Equation 2):
$$L[y] = \frac{d^n y}{dt^n} + p_1(t) \frac{d^{n-1} y}{dt^{n-1}} + \dots + p_{n-1}(t) \frac{dy}{dt} + p_n(t)y = g(t)$$ (2)
Operator Notation and Constant Coefficients
The complexity of $n$ derivatives is consolidated into a single linear operator $L$. When the coefficients are constants ($a_n$), the expression simplifies to:
$L[y] = a_0y^{(n)} + a_1y^{(n-1)} + \dots + a_{n-1}y' + a_ny = g(t)$
This notation emphasizes that $L$ acts linearly: $L[c_1y_1 + c_2y_2] = c_1L[y_1] + c_2L[y_2]$. This principle ensures that the general solution is composed of a complementary solution ($y_c$) and a particular solution ($Y$).
Consider Figure 4.2.4: A two-spring, two-mass system with masses $m_1, m_2$ and displacements $u_1, u_2$. Physics yields two coupled 2nd-order equations. By isolating $u_1$ through substitution, we generate a single 4th-order equation. To solve this, we require 4 initial conditions (position and velocity for each mass) to find a unique physical path.
Worked Example: The Homogeneous Solution
Find the general solution of the differential equation: $y''' - y'' - y' + y = 0$
Assume $y = e^{rt}$. Substituting into the ODE gives: $r^3 - r^2 - r + 1 = 0$.
Factor by grouping: $r^2(r - 1) - 1(r - 1) = 0 \implies (r^2 - 1)(r - 1) = 0$.
This expands to $(r - 1)(r + 1)(r - 1) = (r - 1)^2(r + 1) = 0$.
Roots are $r = 1$ (multiplicity 2) and $r = -1$. Because $r=1$ repeats, we multiply the second term by $t$.
$y_c(t) = c_1e^t + c_2te^t + c_3e^{-t}$